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Abstract 

\ In view of its relation to Big-Bang Nucleo-Synthesis and a reported discrepancy be- 

tween nuclear models and data taken at S-Dalinac, electro-induced deuteron break-up 
£N) . 2 H(e, e'p)n is studied at momentum transfer q < 100 MeV and close to threshold in 

the low-energy nuclear Effective Field Theory without dynamical pions, EFT(^). The 
result at N 2 LO for electric dipole currents and at NLO for magnetic ones converges 
■ order-by-order better than quantitatively predicted and contains no free parameter. 

It is at this order determined by simple, well-known observables. Decomposing the 
triple differential cross-section into the longitudinal-plus-transverse (L+T), transverse- 
transverse (TT) and longitudinal-transverse interference (LT) terms, we find excellent 
agreement with a potential-model calculation by Arenhovel et al using the Bonn po- 
tential. Theory and data also agree well on (Jl+t- There is however no space on the 
theory-side for the discrepancy of up to 30%, 3-<t between theory and experiment in 
(JLT- From universality of EFT(^), we conclude that no theoretical approach with 
the correct deuteron asymptotic wave-function can explain the data. Un-determined 
short-distance contributions that could affect olt enter only at high orders, i.e. at the 
few-percent level. We notice some issues with the kinematics and normalisation of the 
data reported. 
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1 Introduction 



The deuteron is the simplest nucleus, playing the same fundamental role in Nuclear Physics 
as the hydrogen atom in Atomic Physics. Its electromagnetic properties have been studied in 
great detail both theoretically and experimentally - amongst others in photo-disintegration 
(see e.g. [1,2] for reviews) and electro-disintegration (see e.g. [3,4] and references therein). 
The latter process has the advantage to allow for an independent variation of energy and 
momentum transfer. Most experiments have been performed at higher energy transfers and 
are in good agreement with potential-model calculations. However, one experiment [5, 6] 
at the S-Dalinac accelerator of TU Darmstadt (Germany) examined in 2002 the triple- 
differential cross-section for d(e,e'p)n at low momentum transfer (< 60 MeV/c) and close 
to the breakup threshold, concentrating on the decomposition into the different structure 
functions. For the longitudinal-transverse interference cross-section ctlt, a discrepancy of 
about 30% or up to 3 standard-deviations was discovered relative to the prediction by Aren- 
hovel et al. [4,5,7]. As the discrepancy is less pronounced with increasing momentum and 
energy transfer, this does not necessarily contradict other experiments at higher transfers 
which are in good agreement with the same potential-model calculations, see e.g. [8]. 

The disagreement raises however a serious question: The same reaction at the real- 
photon point, i.e. photo-induced low-energy deuteron breakup and recombination np <-» d'y, 
is the first and very sensitive step in Big-Bang Nucleo-Synthesis BBN. Indeed, both the 
higher-energy regime of BBN-relevant photo-dissociation and the S-DALINAC experiment on 
electro-dissociation are most sensitive to the same electric dipole transition amplitude Ely. 
Due to the great difficulty of very-low-energy experiments, the lack of data for this reaction 
in the BBN-relevant energy-region EE [20; 300] keV makes Nuclear Theory at present the 
sole provider of input for the BBN network codes [9,10]. While this has also triggered further 
experiments at S-Dalinac [11] and the High-Intensity 7-ray Facility HI7S at Tunl [12,13], 
the accuracy claimed by theoretical approaches [1, 14-16] is not yet matched. In view of 
this, the S-DALINAC findings are the more troubling: How reliably does theory understand 
the simplest nucleus at low energies, if one of the conceptually simplest non-trivial processes 
seems to disagree with data? Does this reflect deficits in our understanding of the long- 
range NN force, meson-exchange currents, "off-shell effects" and cutoff-dependence, or even 
gauge-invariance? The analysis in Ref. [5] notes that the potential-model approach "without 
meson-exchange, isobar and relativistic currents accounts for 99% of the final result for 
<Jlt, leaving no room for further improvement" and continues: "At present, there exists no 
explanation of this surprising result in the framework of conventional nuclear theory. It is an 
open question whether an alternative interpretation can be offered by effective field theory 
[. . . ]". In this work, we answer this question by concluding that no consistent theoretical 
approach can accommodate the data of Ref. [5]. 

We employ pion-less Effective Field Theory EFT (ft) in the variant in which the effective 
range is re-summed into the two-nucleon propagator [17-23]. EFT(f) is well-tested in low- 
energy reactions with up to three nucleons, see e.g. [24-27] for reviews. Even properties of 
4 He [28-30] and 6 Li [31] are now studied successfully. In particular, np — > was studied 
due to its relevance for BBN [15,32,33], culminating in a N 4 LO calculation by Rupak [14] 
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with a theoretical accuracy of <1%. We here study the S-Dalinac data [5] within the 
same framework, noting that the lowest energy data show the biggest discrepancy and 
lie well within the range of applicability of EFT(/). At momenta below the pion-mass, 
probes cannot resolve the long-range part of the TVTV-potential as nonlocal. Thus, the most 
general Lagrangean compatible with the symmetries of QCD can be built out of contact 
interactions between nucleons as the only dynamical effective hadronic low-energy degrees 
of freedom. EFT(ft) is the unique low-energy version of QCD and allows often (as in this 
case) for simple, closed-form results. The amplitude is expanded in a small, dimension-less 
parameter Q = the ratio between a typical low-momentum scale p typ of the process 
involved, and the breakdown scale A^, set by the mass m n of the pion as the lightest 
particle which is not included as a dynamical effective degree of freedom. This allows one 
to estimate the theoretical uncertainties which are induced by neglecting higher-order terms 
in the momentum expansion of all forces. With a typical low-momentum scale set by the 
inverse deuteron size, one obtains Q « |. Looking at the order- by-order convergence of all 
applications considered so far, one finds this estimate to be quite conservative and the actual 
convergence pattern to be much better. Dimensional regularisation is employed in the two- 
nucleon system to renormalise high-momentum parts of loops, preserving all symmetries at 
each step. External currents and relativistic effects are included in a systematic, manifestly 
gauge-invariant way. No ambiguities arise from meson-exchange currents or off-shell effects. 

The presentation is organised as follows: After recalling the essentials of EFT(^) in Sect. 2, 
the cross-section and individual contributions are given for deuteron electro-disintegration 
up to next-to-next-to-leading order (N 2 LO) for electric dipole transitions Ely and up to 
NLO for magnetic is- vector dipole transitions M\y. Section 3 contains the detailed compar- 
ison with data, focusing on: kinematics in the experiment (Sect. 3.1); the full cross-section 
(Sect. 3.2) and its decomposition into the structure functions (Sect. 3.3); data normalisation 
(Sect. 3.4); and possible higher-order effects (Sect. 3.5). After discussing an experiment at 
slightly higher energies [8] in Sect. 3.6, the Conclusions of Sect. 4 are followed by an Ap- 
pendix with some useful loop integrals and the explicit form of the hadronic matrix elements. 
Further background and details are available in S. Christlmeier's Diplom thesis [34]. 



2 Framework and Calculation 
2.1 Kinematics and Cross-Section 

The kinematic variables of the disintegration process d(e, e'p)n are illustrated in Fig. 1. Two 
frames of reference are conventionally used, see e.g. [3,4,35]. 

On the one hand, the leptonic part e — > 7V of the process is conveniently discussed in 
terms of variables in the laboratory frame (the deuteron rest frame) , denoted by the super- 
script "lab". The four-momenta are (E l y h , k lah ) for incoming, and {E l f°, k' Xab ) for outgoing 
electrons. The scattering angle between the momenta k Ub and k' lab of the incoming and 
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Figure 1: (Colour online) Kinematics of deuteron electro-disintegration, reproduced with 
kind permission of [5]. The electron kinematics refers to the lab frame, while the proton 
variables are defined in the centre-of-mass frame of the two-nucleon final state. 



outgoing electron is Q~ . The energy and momentum transfer of the virtual photon is 

£lab ;= ^lab = £lab _ ^lab ? ^lab = £lab _ £/lab _ ^ ^ 

The hadronic part of the process, 7*d — > pn, is on the other hand more conveniently 
calculated in its own rest frame, i.e. in the centre-of-mass frame of the outgoing nucleons 
(and thus of the virtual photon and deuteron). Variables in this frame carry no superscript. 

2 

The outgoing proton (neutron) has momentum p (— p) and kinetic energy 2M, p '■= \p\- The 
scattering angles in this system are defined by p- q = pq cos Q p and p- (qxk) oc sin $ p . Thus, 
O p is the angle between the virtual photon and the proton, and $ p is the angle between the 
scattering plane, spanned by the incoming and outgoing electron momenta, and the reaction 
plane, spanned by the proton and photon momenta. 

A boost along the momentum transfer q transforms between the two frames: 

^lab J 

/9= M d + u^ ' 7 = 7T^ ' (2 ' 2) 

where Md = 2M — B is the deuteron mass. Hence, the quantities defined in the lab frame 
in (2.1) transform to the proton-neutron cm frame as: 

u = 7 cj lab - /5 7 g lab , q = (3 1 M d . (2.3) 

We choose as the five independent variables of this process the energies E^ h , E l f° (or, 
equivalently, E l ^° = E l ^ h — E l ^ b ) and the scattering angle 0g ab of the electrons in the lab 
frame, plus the proton emission angles 6 P and $ p in the proton-neutron cm frame. The 
momentum of the outgoing proton is 



M 

p= x /(uj-B)M + qi— . (2.4) 
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(0,0,l)g lab , 


Q = 


(0,0,1)9 , 


£lab = 


(sin 6^,0, cos eg')*"' , 


k' lah = 


(sin 9^,0, cos 9 Iab )A; ,lab , 


P = 


(sin 9 P cos $ p , sin 9 P sin $ p , cos 9 p )p 



The momentum vectors are therefore parameterised such that q defines the z-direction, 

(2.5) 
(2.6) 
(2.7) 

(2.8) 
(2.9) 

with the angles between the photon and the incoming or outgoing electrons translating as 

, h k lab - k' Uh cos e!, ab ^ lah A; lab cos 9* ab - A; /lab 

= ^ , cos9> ab = . (2.10) 

Since the transformation between these two frames is given by a boost along q, it does not 
affect the azimuthal angle $ p = <3>„ ab . 

The electron momenta are kf ah = -Eoiab — m 1 an d Mab = -^eiab — m 1- While the electron 
mass m e = 0.511 MeV does not play a role for the experiment at hand, we keep in mind ap- 
plications to corners of phase-space where its effects could be felt, e.g. in back-scattering [11]. 

The triple-differential cross-section is then derived from the amplitude A: 

d 3 a k' lab pM d M 2 



dE l * h dfi' ab dQ p 8 ( 27r )5( M + j£_)J(M d E}? h + <f ab • fc lab ) 2 - M 2 d m 2 e 



\A 2 (2-11) 



d 3 

(TL + &T + &LT cos $ p + (JTT cos 2$ p ) , (2.12) 



d£i ab dfi lab dQ 



where the spherical angles are f2 lab = (9 lab , $|f b = 0) and f2 p = (9 P , <3> p ). In the second 
line, the dependence on the azimuthal proton emission angle $ p is decomposed into the 
longitudinal-plus-transverse parts gl + o~t and the interference terms ctlt and ott [3,4,35]. 
For any process to contribute to either of the interference terms, it must from (2.9) depend 
on those components of the outgoing proton momentum p which are transverse to the photon 
momentum q. This will become important in the analysis later. 

The scattering amplitude A is rewritten in terms of the photon propagator Dj?J, the 
hadronic current J^ adr to be calculated in EFT(^f) below, and the leptonic current / M whose 
contribution to lowest order in the fine-structure constant is easily found: 

A = ^-D^^hadr • (2.13) 

Current conservation implies = = q^J^ adr and is at this order in the fine-structure 
constant a equivalent to gauge invariance. 

If the variables of Ref. [5] are interpreted as discussed in Sect. 3.1, the S-Dalinac 
experiment was performed at two sets of incident electron energies E^f 3 G {50; 85} MeV, 
a range of photon energies E l $° G [8; 16] MeV, electron scattering angle 9^ ab = 40° and 
azimuthal angle $ p = 45°. This leads to photon momentum transfers q G [32; 65] MeV and 
outgoing proton momenta p G [74; 106] MeV in the proton-neutron cm-frame. These are 
the relevant external low-momentum scales of the hadronic matrix element. 
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2.2 Lagrangean and Parameter Fixing 



EFT (ft) allows one to address the questions raised by the S-Dalinac data in a model- 
independent, systematic setting with a reliable estimate of the theoretical uncertainties, 
free of ambiguities, like off-shell effects and cut-off dependence. As the parts of the effective 
Lagrangean of EFT(^), its power-counting rules and the parameter- fixing necessary for this 
work have been discussed extensively by Beane and Savage [22], we repeat them here only 
briefly. The Feynman rules are also given in [34, App. A]. 





= iVt 


iD - 


D 2 
2M 




= -4 


iD 


D 2 
+ AM 




= -A 


iD 


D 2 
+ AM 




C s d 


^ix^jy 




VMpd 



A, 



[kq + KiT 3 )a ■ B 



N 



C 



St 



MyJr Q p d 



A< 

t\s 3 Bi + H.c. 



Sa - y, 
U - Vt 



s{N T P^ So) N + H.c. 
t]N T pf Sl) N + H.c, 



(2.14) 
(2.15) 



t\(N T O x%J N) + }l.c. 



Cq 
Mp, 



-t\[iD ,O tJ ]t 3 ,(2.16) 

i 

(2.17) 



The one-nucleon Lagrangean Cn of the nucleon iso-doublet N = (^) with isospin-averaged 
mass M = 938.9 MeV consists of two parts: First, the kinetic term with minimal substitu- 
tion, = <9 M + ieQA fl , where Q = |(1 + r 3 ) is the nucleon charge matrix and a = |- — j|= 
the fine-structure constant. Second, the interaction via the iso-scalar (iso-vector) magnetic 
moments, hq = 0.44 (k\ = 2.35) in nuclear magnetons. The spin (iso-spin) Pauli matrices 
with vector (iso-vector) index i — 1, 2, 3 (a — 1, 2, 3) are denoted by a-i (r ). 

The Lagrangean for the auxiliary di-baryon field s a (tj) in the 1 S ( 3 Si) channel is C s 
(Ct) [17-19,36]. As we probe only the np system, = + 'ieA ll for both di-baryons. The 

projection operators are Pa = ^gcx2T2T a , Pi ^ = -^oiOi^- The auxiliary-field parame- 
ters are matched to the effective- range expansion of NN scattering such that the effective 
range is re-summed to all powers. This speeds up convergence and simplifies calculations of 
some higher-order effects [17-23]: 



Vs 



A, 



Mr r 



1 



Vt 



■5 7T 



A, 



Mp d 



7 



2 ' 



p (2.18) 



The parameter p encodes the linear cut-off divergence of individual diagrams, calculated in 
the Power Divergence Subtraction scheme (PDS) version of dimensional regularisation [37, 
38]. It is imperative for self-consistency that physical amplitudes are independent of p. This 
is indeed the case. 

For the spin-singlet, the scattering length and effective range of the np system are probed 
in deuteron electro-disintegration: a = —23.71 fm, r = 2.73 fm. Parameters for the spin- 
triplet are determined from the effective- range expansion about the observed real bound 
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state, namely the deuteron with binding energy B = 2.225 MeV (corresponding to a mo- 
mentum 7 = V MB = 45.70 MeV) and effective range p d = 1.764 fm. Therefore, the 
deuteron wave function shows the correct exponential decay and normalisation Z already 
at LO in EFTfof) [22,39]: 



^deuteronO ~> Oo) 



Z e 



with Z 



IPd 



The last two terms of C t parameterise SD-wave mixing in the spin-triplet, with 



O 



1 



%y,3 



l [D X D y P J 



( 3 Si) 



p 



( 3 Si) 



D X D V - D x P, 



( 3 5i) 



D, 



DyP, 



( 3 Si) 



O 



- I DiDj 



3 3 



(2.19) 

(2.20) 
(2.21) 



operators which ensure that C s t is manifestly gauge-invariant. Their strengths are deter- 
mined from the asymptotic ratio rj sd = 0.0254 of D- and S-wave components of the deuteron 
wave function, and of the deuteron quadrupole moment \iq = 0.2859 fm 2 [22]: 



G 



sd 



Mr 



2Z 



Vt 



C sd M 2 2 



C c 



v ^7^327r3 7 M Pd \ 



(2.22) 



Here, C sd contributes at LO, and Cq provides a NLO correction of about 50%. 

Finally, C st in (2.17) parameterises transitions between the 1 Sq and 3 Si channels by 
a magnetic field acting on a di-baryon. Its strength is determined by the thermal cross- 
section <j(E = £ = 1.264 • 10~ 8 MeV) = (334.2 ± 0.5) mb [40] for radiative capture of 
neutrons on protons np — > dj. At thermal energies, it is dominated by iso-vectorial M\y 
transitions. Magnetic-moment interactions are LO, and the parameter L\ in £ st enters as a 
NLO correction of about 50%, while electric transitions are irrelevant. The free parameter 
is thus determined from the thermal cross-section [22] 



(7 



(Afl v ) 



2aZ(p 2 + 7 
pM 3 



2\3 



1 



1 9 
2 r 0P 



ip 



7 ~ \ + \roP 2 U 
1 p 2 + 7 2 + 2 



(2.23) 



as Li = —4.41 fm [34] or L\ = —4.0 fm [22], depending whether terms quadratic in L\ are 
kept. Both variants differ by 10% < Q, as expected from L x being a higher-order correction. 

We now state the power-counting rules of EFT(^) in the version in which the effective- 
range parameters are treated as unnaturally large and thus are kept at LO together with 
the scattering lengths [22]. The typical low-energy centre-of-mass momentum p and kinetic 
energy E/2 of a nucleon is counted as 



\p\ ~ Vme 



1 1 1 

Pd r a 



p ~ QAf 



(2.24) 



Hence, the wave function renormalisation scales as Z ~ 1. SD mixing is suppressed by Q 2 
with respect to pure S'-wave amplitudes because of the two derivatives in (2.16). Finally, 
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the ratio of the iso-scalar and iso- vector magnetic moments is ^o/ K i ~ 1/5 ^ Q ■ Therefore, 
fto is neglected as numerically higher order in the following [41]. 

Since the gauge field is minimally coupled, these rules are simple to extend: The 
zero-component of the gauge field scales like an energy, A ~ Q 2 ; its 3- vector com- 
ponents like a momentum, A ~ Q. As seen in the previous Sub-Section, photon en- 
ergies for the S-Dalinac kinematics of deuteron electro-disintegration lie in the range 
E l * h G [8; 16] MeV, corresponding to E x G [6; 14] MeV ~ [Q 2 A^/M; Q 2 A^], photon mo- 
menta in the range q G [32; 65] MeV ~ QA^, and outgoing proton momenta in the range 
p G [74; 106] MeV ~ QA^, approaching A^ « m n . Convergence must therefore be checked 
with care. 



2.3 Electric Contributions to N 2 LO 

Electric dipole transition amplitudes Ely dominate not only the photo-dissociation cross- 
section at the higher end of the energy range relevant for BBN, but also electro-disintegration 
for p > 20 MeV. To N 2 LO, the hadronic currents have the structure 

= ^^(K^K^Jtj , (2.25) 

where the proton and neutron spinors N v / n are specified explicitly. The indices on the 
current indicate the polarisation j of the initial deuteron and the spin i of the pra-triplet 
final state. Averaging over spins and polarisations, the electric contribution to (2.11) is 

\A Elv \ 2 = ( { T a) \s£ [(*& Jj* ^b4 j + (A* - ")) - (^ab - ml)J^J^) . (2.26) 
{u — q ) o 

The three diagrams contributing at leading order, Q°, are shown in Fig. 2. In the 




Figure 2: (Colour online) The LO electric contributions to the hadronic current. The double 
line denotes the spin-triplet di-baryon intermediate state. There is no NLO contribution. 



last two graphs, the photon couples to the deuteron which is at this order a pure S-wave 
state. These diagrams are therefore independent of the direction of the outgoing proton 
momentum p, and hence via (2.9) of the proton emission angle $ p . They do thus not 
contribute to the longitudinal-transverse and transverse-transverse parts of the cross-section 
(2.12). In contradistinction, the final-state interaction in diagram (el) depends on p and 
thus contributes to olt and ott- 
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(sdl) (sd2) (sd3) (sd4) 




(sd5) (sd6) (sd7) (sd8) 



Figure 3: (Colour online) The order-Q 2 (N 2 LO) electric contributions. The SD-mixing 
vertex proportional to Csn is denoted by a circle; further notation as in Fig. 2. 



There are no electric contributions at NLO. The first non-zero corrections, listed in Fig. 3, 
appear at N 2 LO. All arise from SD-mixing proportional to Csd in the Lagrangean (2.17). 
Relative to the LO contributions, these diagrams are suppressed by two more derivatives, 
i.e. two more powers of Q « ~ and thus should contribute at most (~) ~ 10% of the LO 

terms. The error one makes by truncating the series at this order should be (|) 3 « 4%. 
The convergence is indeed much better because the asymptotic SD-ration is numerically 
rjsd = 0.0254 ~ Q 3 , i.e. N 3 LO. Even at the S-DALlNAC-kinematics of photon momenta as 
large as 78 MeV (E^ = 50 MeV, E l f = 9 MeV, 6^ b = 40°, $ p = 45°), they amount to 
not more than 1% of the total; see the discussion and figures in Sects. 3.2 and 3.3 as well 
as Ref. [34] for details. Only the diagrams (sdl-3) and (sd6) contribute to the interference 
terms <jlt and &tt via their final-state interaction or via the -D-wave component induced 
into the deuteron and np final-state wave-function. 

The resulting electric currents up to N 2 LO are detailed in App. A. 2. At each order, they 
are manifestly gauge-invariant and cutoff-independent. 



2.4 Magnetic Contributions to NLO 

Magnetic dipole transitions Mly dominate the inverse reaction np — > d'j at the lower end 
of the energy range relevant for BBN, but are usually small at S-Dalinac kinematics. We 
therefore include them only to NLO, Fig. 4. The LO part is set by the iso-vector magnetic 
moment K\\ the only NLO graph, (m4), comes from the singlet-triplet coupling (2.17). 

The hadronic currents from magnetic contributions can to this order be written as 

& " = eVz^^Ni^K) J 3 , (2-27) 
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Figure 4: (Colour online) The magnetic contributions to the hadronic current at LO (ml- 
m3) and NLO (m4). The blob denotes a magnetic photon coupling via the iso- vectorial 
magnetic moment; the square the dibaryon-photon interaction (2.17) proportional to L\\ 
the dashed double-line the spin-singlet di-baryon intermediate state; further notation as in 
Fig. 2. 

and the squared amplitude is after averaging over initial states: 

\Ami v ? = A 47ra L f (S jm S kn - 5 km V n ) (2.28) 
(or — q z ) z 6 

x [k'^ h J ]j kf h J m + k'^ h J m k^ h J ]j + (k^ h k' 1 ^ - ml) J tj ; J m S kn ] 

The vertex where a magnetic photon couples to a di-baryon reduces the amplitude largely, 
although it is formally NLO. Compared to the electric amplitude, the magnetic one is 
suppressed by two orders of magnitude in the energy regime of the S-Dalinac experiment, 
except for the vicinity of Q p = 70°; see Figs. 5 to 8. 

It is not too difficult to see that the only angular dependence in the magnetic amplitude 
comes from the p ■ q -terms in the intermediate nucleon propagators, as the deuteron and 
outgoing np-system are again pure S-states in these diagrams. Thus, the amplitudes only test 
the longitudinal part of p. There is no $ p -dependence, and therefore magnetic transitions 
do at this order not contribute to o^t and ott interference terms (2.12). 

The explicit form of the resulting magnetic currents can be found in App. A. 3. At each 
order, they are manifestly gauge-invariant and cutoff-independent. 

3 Theory Confronts Data 

We now compare the results of EFT(f) to the S-Dalinac data reported in [5,6] and to the 
potential model calculation by H. Arenhovel et al. which is based on the Bonn-C r-space 
potential and includes final-state effects, meson exchange currents, isobar configurations 
and (in our case negligible) relativistic effects [4,5,7]. Before we consider the energy- and 
angular dependence of the triple differential cross-section and its decomposition, we first 
have to address a subtle kinematical point. 

3.1 Kinematics, Again 

In Sect. 2.1, we adopted the standard kinematics of electro-disintegration according to which 
the leptonic part of the cross-section is calculated in the lab frame (deuteron at rest), while 
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the hadronic part is determined in the proton-neutron cm-frame of the hadronic sub-process 
7*d — * pro. In a slight but revealing abuse of language, the latter is often referred to as "the 
cm frame". According to a literal reading of Ref. [5], the S-Dalinac experiment would have 
been analysed in the centre-of-mass frame of the whole process, resulting in large deviations 
between theory and experiment both in shape and normalisation even of the full triple- 
differential cross-section, see Fig. 5. We take this to be a slip of the tongue and assume that 
the experiment was analysed using standard kinematics [3,4,35]. 

Furthermore, the variable which labels the energy transferred by the photon in the S- 
Dalinac experiment is said to be obtained "after transformation to the centre-of-mass 
system" [5]. However, agreement of the data with both the potential-model and EFT(^) 
calculation is found only if the S-Dalinac- variable "Ex" refers to the lab frame. The point 
£jPn)-cm = 9 MeV would correspond to E 1 ^ = 10.6 MeV and a much smaller cross-section, 
so that the resulting curve would again substantially deviate from data, see Fig. 5. We 
therefore interpret "Ex" in [5] as E^, defined in the rest frame of the deuteron target. 

The experimental results are reported in 4 energy bins, E 1 ^ = [8; 10], [10; 12], [12; 14] and 
[14; 16] MeV. We interpret them such that the count-rate was normalised to give the double- 
differential cross-section per MeV. Thus, we compare the triple-differential cross-section for 
E l x h to the double-differential one per MeV. The difference between averaging over the 
2 MeV-range and taking the mean value of E 1 ^ can be neglected down to = 6 MeV, as 
shown in the plot; see also [34]. If Ex is interpreted as given in the pn-cm frame and the 
cross-section as given over the whole 2 MeV-bin, the data would consistently over-shoot the 
theory results by 10 to 20%. 

We are grateful to H. Arenhovel [7] for confirming that the potential-model results re- 
ported in [5] were determined in the same kinematics used here. 

3.2 Total Differential Cross-Section 

Comparing in Fig. 6 to the differential cross-sections reported in Ref. [5] at E}f° = 50 MeV 
and 85 MeV in several i±7^ b -bins confirms our interpretation of the experiment's kinemat- 
ical variables. At E}^ h = 85 MeV, additional data at large Q p are available and within 
their (combined statistical and systematical) error-bars compatible with both theoretical 
approaches. We note excellent agreement with the calculation by Arenhovel et al. and with 
the data within error-bars. In Fig. 6, a decomposition of the different contributions in 
EFT(jf) for the smallest and the largest energy and momentum transfer of the S-Dalinac 
data is also shown. As expected, the LO electric transitions dominate except around the 
minimum at p ~ 70°, where NNLO contributions and magnetic transitions play a sig- 
nificant role. That means that by taking into account only minimal coupling of photons 
to nucleon and di-baryon fields, disintegration can be described highly accurately within 
EFT (ft). SD mixing, being the only correction up to NNLO, contributes even less than es- 
timated by power counting. This good convergence confirms our confidence that our results 
are reliable. Less-known short-distance contributions play a very minor role, see Sect. 3.5. 

The convergence plots also reveal an interesting facet: The N 2 LO contribution has in- 
creased slightly for higher energy and momentum transfer, but much less than expected 
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Figure 5: (Colour online) Examples of the triple-differential cross-section for E l £ h = 9 MeV 
at E^ h = 85 MeV, compared to S-Dalinac data [5] (with combined statistical and system- 
atical error-bars) and to the result by Arenhovel et al. (squares) [4,5,7]. (a): Interpretation 
of "Ex = 9 MeV" as given in the total cm frame, the pn-cm frame, or the lab-frame; (b): 
EFT(/) at E l f = 9 MeV (solid line) and £^ pn) " cm = 9 MeV (i.e. Ef> = 10.6 MeV, dashed). 
The dotted curve is the double-differential cross-section per MeV, obtained by integrating 
over the bin G [8; 10] MeV; the difference to the solid line is almost invisible, (c): 
comparison on a logarithmic scale of the EFT(^) contributions: electric transitions up to 
LO and N 2 LO respectively (recall that NLO-ji| zero), and magnetic transitions up to NLO. 



E a = 50MeV,E a = 9 MeV, 9 a =40°, <b =45° E„ a = 85 MeV, E a = 15 MeV, 9 a =40 o ,<b =45° 

U X e P X e P 




Figure 6: (Colour online) EFT (ft) results for the differential cross-section, compared to data 
and calculations by Arenhovel et al. (a): data set at smallest energy and momentum transfer, 
Eq* = 50 MeV, E^ h = 9 MeV. (b): data sets at largest energy and momentum transfer, 
£' ab = 85 MeV, Ef 3 = 15 MeV. (e,f): data set at ^ ab = 85 MeV, Ef> = {11; 13} MeV. 
Comparisons of EFTfjf) contributions are in Figs. (c,d) included for the smallest and largest 
energy and momentum transfer: electric transitions up to LO and N 2 LO respectively (NLO 
is zero), and magnetic transitions up to NLO. 
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from naive dimensional analysis. Deviations between the potential-model and EFT (ft) re- 
sults are also slightly more pronounced there, especially at large 9 P . As typical momenta 
in the process approach the breakdown-scale of EFT (ft), the dimension-less parameter 
Q = Ptyp/Ajf approaches unity and the expansion is rendered useless. The momentum trans- 
fer q G [32; 56] MeV is comparable to the intrinsic low-momentum scale of the two-nucleon 
system, namely the deuteron "binding momentum" 7 ~ 45 MeV. The momentum of the 
final-state proton in the proton-neutron cm frame ranges however from 75 to 106 MeV and 
hence becomes at the upper end comparable to « m n . One would have expected this to 
affect in particular the final-state interaction diagrams and the SD mixing contributions at 
N 2 LO, which depend on p in the numerator. The observed enhancement is however much 
smaller. The good convergence of EFT (ft)- even for momenta close to the breakdown scale 

- was found also in many other applications, see e.g. references in [25-27,42]. 

3.3 Comparing Structure Functions 

The main goal of the S-Dalinac experiment [5] was a decomposition of the contributions 
of different structure functions in (2.12): the longitudinal-plus-transverse cross-section, and 
the LT- and TT-interference parts. However, the uncertainties in the data turned out to 
be too large for a meaningful comparison of ott to theory. Studying the interference terms 
provides information about the impact of final-state interactions, since these can be "In- 
dependent, as re-iterated in Sect. 2.1. The contributions of the L + T, LT and TT terms are 
shown in Figs. 7 and 8 for E^f° = 85 MeV and E l ^° = {9; 11; 13; 15} MeV, again compared 
to data and the potential-model result. 

The results are normalised to d 3 (<r L + a T ) / (dE l e ah dQ l ^ b dQ p ) at Q p = 0° to make the 
comparison independent of the absolute normalisation of the data. This also mutes the 
question whether the cross-sections were averaged or summed over each _E^ b -bin. The error 
of the value to which the data are normalised has been taken into account in all error bars. 
o\l + °t dominates, whereas and o?t are one and two orders of magnitude smaller, 
respectively. However, the interference terms become more relevant with increasing energy 
transfer, reflecting the growing impact of final-state interactions. In each figure, the result 
for the LO electric transitions in EFT(/) is also given. Recall that in that case, glt and ott 
are nonzero only because of the final-state interaction of the proton with the virtual photon 
after the deuteron breakup, Fig. 2. (el). The degree to which the -Ely-transition dominates 

- already discussed in the previous Section - is even more striking for ctlt- Magnetic 
transitions do not contribute at all to the interference cross-sections, since their amplitudes 
are independent of $ p , see Sect. 2.4. The impact of SD mixing is increasing slightly with 
E l x h , but is small for a L + a T and <ttt, and almost negligible for ctlt- 

The EFT (ft) results are again in good agreement with those of Arenhovel et al. Especially 
Olt is essentially identical in both calculations. Only in the very small quantity Ott can a 
deviation be detected. It is a factor of 4 larger than the difference between the LO and N 2 LO 
results, which in turn is one way to estimate higher-order effects. However, we already hinted 
above that the contributions from SD mixing and magnetic moment interactions are much 
smaller than estimated by naive power-counting. Another estimate for possible corrections is 
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Figure 7: (Colour online) Decomposition of the triple-differential cross-section into (top to 
bottom) longitudinal-plus-transverse (L + T), longitudinal-transverse (LT) and transverse- 
transverse (TT) parts at E^f 3 = 85 MeV, normalised to (Tl+t at P = 0. Left column: 
Ex h = 9 MeV; right column: E l ^° = 11 MeV. No data are available for the TT interference 
cross-section. 
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E lab = 85 MeV, E '* = 13 MeV, 6 bb = 40° 



E„ lab = 85 MeV, E = 15 MeV, 6 lab = 40° 
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Figure 8: (Colour online) As Figure 7, for E% b = 13 MeV (left column); E% b = 15 MeV 
(right column). 
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to set the uncertainty of the EFT(^) calculation at N 2 LO as ~ Q 3 = (p typ /A^) 3 « [3 . . . 10]% 
of the LO term, assuming a typical momentum p typ ~ q ~ [45; 60] MeV. This is apt to 
over-estimate the theoretical uncertainties, as discussed in Sect. 3.5 in conjunction with 
possible higher-order corrections. In short, all EFT(^) results should be understood with 
an uncertainty which reflects an estimate of the higher-order effects not considered. A 
conservative accuracy limit would be < 10% in each observable. It is thus safe to conclude 
that the potential-model and EFT(ff) results agree within the theoretical uncertainties. 



3.4 How To Resolve Discrepancies with Data? 

The accordance of both theoretical approaches with the data, however, is considerably 
worse. The discrepancy ranges from just over one standard deviation at the largest energy 
and momentum transfer, to three and more at the smallest energy and momentum transfer. 
Not only does the calculated longitudinal-transverse interference cross-section deviate by 
about 30% from the measured one (as pointed out in [5]), but also a L + o T at large 9 P . In 
particular, it is remarkable that the latter is over-predicted at E l ^° = 9 MeV, but under- 
predicted at E l ^° = 15 MeV. A look at <jl+o~t at 6 P = 180° before data normalisation reveals 
however a reduced discrepancy, see Fig. 9. Only then do data and theory agree on the shape 
of the curve, except at the smallest value of E l £ b , and an additional normalisation would 
solve the problem. While the statistical significance of the discrepancy hardly changes, this 
may indicate un-accounted systematic errors in the data. 

Normalising to P = seems to deteriorate the accord of calculations and data - although 
it avoids accounting for the absolute normalisation of data. The reason can be seen by 
comparing Fig. 9.(b) and (c): The cross-sections both at P = and P = 180° display 
relatively small deviations which go, however, in different directions: The theoretical result 
e. g. at Ex h = 15 MeV is a bit too small for O p = 180°, but too large for P = 0. In dividing 
these two values, the discrepancy increases. 

The significant difference between data on the one hand, and the agreeing theoretical pre- 
dictions of EFT(ff) and Arenhovel et al. on the other hand, continues for the longitudinal- 
transverse interference cross section glt, see Figs. 7.(b,e) and 8.(b,e). The angular depen- 
dence of the normalised values is reproduced, but the normalisation differs between one 
standard deviations at high energy and momentum transfer, and three standard deviations 
at high energy and momentum transfer. 

Figure 10 examines the energy dependence of the discrepancy for o^ T at the minimum 
around p = 35° as function of E^ h , both before and after normalisation of the data. 
Again, the situation becomes, at least for higher energies, slightly better when looking at 
the absolute values. However, a clear deviation remains, and surprisingly becomes larger at 
smaller E l ^ h . The shapes and slopes of the un-normalised energy-dependence of the minima 
in ctlt do not match. The un-normalised data at high energy are, at least at non-backward 
angles, compatible with the theoretical approaches within error-bars, but the discrepancy 
at smaller energies and momentum transfers worsens, see Fig. 11. 
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e = 180° 




e =o 




Figure 9: (Colour online) £jf b -dependence of ol + &t at P = 180°, normalised to Q p = 0° 
(a) and before normalisation (b). Bottom: + ut at 6 P 
normalised). All graphs refer to E l Q ab = 85 MeV, Of = 40°. 



0° (to which all data are 



3.5 Discussion of Higher-Order Contributions 

In view of these discrepancies, one must investigate to which extend additional interactions 
could remedy the problem. The EFT(^) result is complete up to N 2 LO in electric and up to 
NLO in magnetic transitions, resulting in a conservative accuracy-estimate of < 10%, as dis- 
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Figure 10: (Colour online) E£ -dependence of o~lt at the minimum around Q p = 35°, 
normalised (left) and absolute values (right), both for E l ^° = 85 MeV, 0^ h = 40°. 
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Figure 11: (Colour online) Ej£ -dependence of the absolute values of ctlt for the extreme 
cases E l f = 9 MeV (left) and 15 MeV (right) at E l ah = 85 MeV, 9^ b = 40°. Compare to 
the normalised data in Figs. 7 and 8. 



cussed in Sect. 3.3. Including some renormalisation-group invariant sub-set of higher-order 
contributions provides a third alternative to estimate the uncertainties of the calculation, 
besides order-by-order convergence and prima facie power counting. 

However, the results so far provide a strong argument that higher-order corrections 
cannot solve the problem: Recall first that a^x is nonzero only because of non-zero transverse 
components of the momentum p of the outgoing proton. If, in spite of the power-counting, 
some higher order contributions were important for (Jlt, they should therefore contribute 
more at higher energies, where p is bigger and thus closer to the break-down scale A^. 
This would lead to larger higher-order corrections with increasing photon energy. This is 
however at odds with our analysis of Figs. 7 to 10. As fine-tuning could circumvent this 
general argument, we explicitly consider in the following some N 3 LO contributions. 
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Magnetic transitions are negligible at these energies and in these kinematics, as seen in 
Figs. 5 and 6. Another contribution to electric transitions includes P-wave nucleon-nucleon 
final state interactions which also enter at N 3 LO [33]. While they do become stronger with 
increasing energy and momentum transfer, they are, even at high energies, not even of the 
same order of magnitude as the N 2 LO contributions. 

Next, one might think that the correction to the deuteron quadrupole moment of about 
50% provided by the Cg-term in (2.16) can solve the problem at N 3 LO. However, this term 
is at best comparable in size to the N 2 LO term from SD mixing. Its contribution will thus 
be even smaller than the N 2 LO correction in Figs. 7 and 8. An explicit calculation [34] 
renders 

4 = _ 2 y,A^(g^-^q%) , 4 = . (3.1) 

The difference to the N 2 LO result is less than 0.5 % even at the minimum around P = 70° 
at E l Q b = 50 MeV, E l ^° = 9 MeV, where the effect would need to be biggest to remedy the 
discrepancy. The contribution of this interaction to the interference terms ojjf and ott is 
furthermore zero since (3.1) is independent of $ p . 

Corrections from relativistic kinematics are negligible. An example for a relativistic effect 
which is dynamical but has been neglected in our calculation is the spin-orbit interaction 
first considered in EFT(^f) by Chen et al. [43,44], 

C so = iAt 

where E is the electric field. This term is suppressed by p typ /M < Q 2 relative to the 
magnetic-field term in Cn (2.14) and thus enters in M\y transitions at N 4 LO. Chean et 
al. also demonstrated that it provides substantial angle-dependent contributions to deuteron 
Compton scattering at energies around u rs 49 MeV [44]. Its Feynman diagrams are those 
of the LO magnetic transitions, Fig. 4.(ml-3), with the magnetic-moment interaction sub- 
stituted by (3.2). The vector-component of the hadronic current has the same structure as 
(2.27), and the zero-component can be written as 

J°S lv) = e^h\ d) ^{Nla 2 N:)J% . (3.3) 

The analytical results for these currents are reported in [34]. As expected, they give only a 
small correction to the already negligible NLO result for magnetic transitions. In (Jl+t, they 
lead to a correction never exceeding 10~ 3 of the LO result. Although they are $ p -dependent, 
the relative difference to LO in cjlt is less than 10~ 2 , with 2 x 10~ 3 at the extrema. In ott, 
the relative difference to LO at the maximum is 10 -5 and never exceeds 10~ 4 . 

We mention in passing also a N 4 LO contribution which can naively be incorporated by 
modifying the di-baryon propagators (A.12/A.25) to include the shape parameters p\ = 
0.389 fm 3 for the 3 Si channel and r x = —0.48 fm 3 for the 1 S channel of iVTV-scattering. 
This neglects diagrams which come from gauging the derivative terms in the Lagrangean 
which correspond to these corrections, so the result is not complete. However, we see that 
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Figure 12: (Colour online) Comparison to the data of Ref. [8] for o~l+t (a) and o^t (b), 
both in the lab frame (also for the hadronic variables!) and normalised to cr L+T at P = 0. 



this leads only to minimal modifications of the structure functions. (Jl+t is increased by 
< 5% at back-angles, about 20% in the minimum due to fine-tuning, and in general much 
less, ojjt changes by less than 1% at the extrema, and o~tt only by 0.1%. 

Overall, this discussion supports our estimate of the theoretical uncertainty of our results. 

3.6 A Higher-Energy Experiment 

In Fig. 12, we finally compare the EFT{f) results to an earlier experiment at a slightly 
higher momentum transfer [8] which reported good accord between data and the potential 
model calculation by Arenhovel et al. 1 also for the normalised <jlt- Although the involved 
proton momentum p = 117 MeV is even closer to the breakdown scale of EFT (ft) than for 
the S-Dalinac kinematics, we observe very good agreement for both (Jl+t and a LT . It is 
remarkable that the result describes the data so well in a regime in which EFT (ft) may start 
to become unreliable, but deviates significantly from data at lower momenta. 

4 Concluding Questions 

In a parameter-free N 2 LO calculation using EFT(/), a manifestly gauge and renormalisation- 
group invariant, well-tested formulation of few-nucleon physics with analytical results which 
are rooted in a systematic, model-independent low-energy expansion of all nuclear forces, 
with an a priori estimate of theoretical uncertainties corroborated by order-by-order con- 
vergence, the S-Dalinac data on deuteron electro-disintegration d(e, e'p)n at low energy 
and momentum transfer [5,6] cannot be explained. 

1 For this experiment, the calculations were performed using the Paris potential; the difference to the 
Bonn potential calculation in Ref. [4,5,7] should not be relevant here. 
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Since an important aspect of any Effective Field Theory is its universality and result- 
ing model-independence, this leads us to a conclusion beyond the EFT(^) approach: The 
predictions of models which reproduce or share the input of EFT(?f) must agree with its 
results to within the accuracy of the EFT calculation, in the range where EFT(^) is ap- 
plicable. The analysis in Sect. 3 showed that the S-Dalinac data at low momentum and 
energy transfer lie undoubtedly within this regime. Therefore, any self-consistent potential 
model, irrespective of the detail led treatment of meson-exchange currents, "off-shell effects", 
cutoff-dependence, etc. will to the accuracy outlined here agree with the EFT(?f) result, if it 
incorporates the same ingredients: the same deuteron binding energy, the same asymptotic 
normalisation Z and asymptotic S-to-D wave ratio r] sd of the deuteron wave-function, the 
same scattering length and effective range in the 1 S channel, the same iso-vectorial mag- 
netic moment of the nucleon K\ and the same total radiative thermal capture cross-section 
of neutrons on protons. This is nicely confirmed for electro-disintegration at low energies: 
The potential-model approach by Arenhovel et al. [4,5,7] yields essentially the same results 
as EFT(^), especially for glt- 

We find furthermore that cr L+r , <jlt and o^t are all dominated up to the few-percent 
level by the leading-order electric transition, i.e. by minimally coupling the photon to the 
nucleon and deuteron. This coupling is sensitive only to the asymptotics of the deuteron 
wave- function. The statement is strongest for Olt-, which shows the discrepancy: It is 
dominated up to the 1%-level by only one LO electric process, Fig. 2. (el), i.e. by the 
asymptotic normalisation Z of the deuteron wave- function, eq. (2.19). In contradistinction, 
the discrepancy to the S-Dalinac data amounts to up to 3 standard-deviations or 30%. 

These findings suggest a re-analysis of the experiment [45]. We identified some questions 
concerning the kinematics and systematics of the experimental analysis and caution that 
the differences could in part arise from a dis-advantageous normalisation of the data. If 
the discrepancies were confirmed, this would pose a highly non-trivial problem for Nuclear 
Theory. Only new data for deuteron electro-disintegration near threshold and at low mo- 
mentum transfers can settle the issue definitively, particularly concerning the decomposition 
into the contributions of different structure functions. We are however confident to maintain 
that the validity and error-estimate of the .Ely-part of the photo-dissociation cross-section 
relevant for Big-Bang Nucleo-Synthesis as calculated in EFT(^) [14, 15, 32, 33] are not in 
question. 
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A Currents and Integrals 



We now list the analytical results for the hadronic currents described in Sect. 2.3 and 2.4. 



A.l Useful Loop Integrals 

The following loop integrals are calculated by the methodology outlined in [46, 47] 2 , using 
contour integration for the energy part and the PDS subtraction scheme combined with 
dimensional integration in D spatial dimension [37,38] to identify divergences, parameterised 
by /x, and tensorial reduction. The reader may consult [34, App. D] for details. 
The fundamental bubble-sum integral is [37, 38] 




(A.l) 



Reference [46] provides 

t(2), , v f d 3 / 1 1 1 / q \ /A , 

Ik (a, o; a) = / -, — r^-= =; — = = arctan =- . (A. 2) 

From that, one finds the /i-independent result with one loop momentum in the numerator 

A(a,b;q) = J ^ )D ± 'J- (f + ^ - - = if (a) - if (b) - (£ + b - a) if (a, b; q) (A.3) 

and the following /i-dependent ones with two powers of loop momenta in the numerator: 
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2 We use a slightly different notation. 
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Three powers of loop momenta in the numerator are covered by the integrals 
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A. 2 Electric Currents up to N 2 LO 

The electric contributions to the hadronic current at LO are labelled as in Fig. 2. The 
4- vector contribution is for the first diagram with r :— 2f> — q 
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where the proton propagator 
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(A.10) 



could be approximated for real photons by dropping the g 2 -term in the denominator (u = q), 
but not in electro-disintegration because the kinematics imposes uj ~ q 2 /M. 
The second diagram gives 
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with the spin-triplet di-baryon propagator 
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After integrating over the loop, diagram (3) results in 
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Here, as in the following, — p 2 is understood to be — p 2 — ie with e \ 0, i. e. \f—p 2 = — ip, 
where p is the proton momentum. 

The N 2 LO contributions come from SD-mixing operators only. Following the numeration 
of Fig. 3, one finds: 
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The last three diagrams contribute only to the three- vector component of the current: 
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Dt{p)y 2 t ^[l£\l 2 ) - Ii>\-P 2 )} (qAk + QiSjk - \q k 5i^j (A.22) 

With these results, it is simple to show that the LO and N 2 LO electric currents are separately 
gauge-invariant, q^J^ = 0. It is also noteworthy that the total current at each order is 



independent of the regularisation parameter fx. 
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A. 3 Magnetic Currents up to NLO 

The LO diagrams, numbered as in Fig. 4, contribute 

J* ml > = -y t D p (p; u, q)^q , = -y t D n (p; u, q)^ q (A.23) 

j(m,3) = ^D^y^MI^i-p 2 , 7 2 ; g) <f . (A.24) 
with neutron propagator D n (p;u,q) = D p (—p;cj,q) and spin-singlet di-baryon propagator 

iDs(p) = ^ ' . • (A.25) 

2 iS-p 2 + in 

ao 2 r f 

The current provided by the NLO contribution is 

J>4) = _ 2ysjDs(p) _^ - (A 2g) 

Notice that these currents are all transversal only. 
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